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1. Introduction 

Directed polymer in a random medium is a widely studied model in the statistical 
mechanics of disordered systems [TJ [2]. The polymer chain is immersed in a static 
random potential. In the directed version there is a singled out direction, also referred to 
as "time" , such that the polymer is constrained to move forward along the time direction. 
Prominent realizations are domain walls in two dimensional disordered magnets [3] , tear 
lines for paper sheets jl] and vortex lines in disordered superconductors [5] . For spatial 
dimension d < 2, the polymer is super diffusive at any coupling strength, while for d > 2, 
there is a weak coupling phase with diffusive behavior and a still little explored strong 
coupling phase, see [(J [7] and references therein to earlier work. 

Recently there has been considerable progress [8j [H [10l [TlJ [121 [131 E] f° r a 
particular continuum version of the directed polymer in 1 + 1 dimension with both 
endpoints fixed (point-to-point directed polymer). The free polymer is modeled by a 
continuum Brownian motion with bending coefficient 7. Then, in the presence of a 
disorder potential V, the point-to-point partition function of the polymer at inverse 
temperature (3 reads 

Z(x,t)= / V[x(r)]exp(-/3 dr [I 7 («9 T x(r)) 2 + V(x(t), t)]) . (1.1) 

Jx(0)=0 V JO J 

The partition function is the sum over all possible paths x(t) of the polymer, starting 
at position at time and ending at position x at time t. The energy of the polymer 
is the sum of the elastic bending energy, proportional to 7, and the potential energy 
obtained from summing the external potential V along the polymer chain. The partition 
function Z(x,t) is random as inherited from the randomness of the potential V, which 
is assumed to have a Gaussian distribution with mean and covariance 

(V(x, t)V(x', t')> =DS{x- x')5(t - t') . (1.2) 

The particular choice of the covariance fll.2j) allows to express the n-th moment of Z as 
a propagator matrix element of an n-particle attractive £-Bose gas on the line, a model 
which can be solved exactly by the Bethe ansatz. As will be explained in more detail 
below, the progress alluded to refers to an exact computation of the generating function 
for the partition sum Z(x,t). In our contribution, this result will be extended to the 
generating function jointly of Z(xi,t) and Z(x2, t), i.e. for two distinct positions X\ and 
£2 of the endpoint of the polymer at the same time t. 

An additional interest in the continuum directed polymer comes from the connection 
to the Kardar-Parisi-Zhang (KPZ) equation [15J, which is a stochastic evolution for a 
growing surface. If we denote the height profile by h(x,t), then, in the conventional 
units, the one-dimensional KPZ equation reads 

d t h(x,t) = ±\(d x h(x,t)) 2 + udlh(x,t) +ri(x,t) . (1.3) 

Here, A is the strength of the nonlinear growth velocity, v is the parameter governing 
the surface relaxation, and 77 is a white noise with strength ypD modeling the random 
nucleation and deposition events at the surface. 
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The partition function of the directed polymer Z(x,t) given in Eq. (II. ip satisfies 

d t Z(x, t) = ^- d 2 x Z(x, t)-/3 V(x, t)Z(x, t) , (1.4) 

from which it follows that the free energy, defined by 

F(x,t) = —log Z(x,t) , (1.5) 

is a solution of the KPZ equation (11. 3p under the identification 

h = -F, A = 7 -\ u=(2f3 1 y 1 , V = -V. (1.6) 

In the context of surface growth, the joint distribution of the free energy is a natural 
quantity of interest: for given time t, it is the joint height statistics at the two spatial 
reference points x\ and X2, which in particular determines the height-height correlations 
at time t. 

Our considerations are somewhat formal, since taken literally (Z(x,t)) = oo. In 
dimension d — 1, this divergence can be easily taken care of by a suitable free energy 
renormalization, as discussed in detail in |16|. After renormalization one finds that 



(Z(x,t)) = (/3 7 /2^r 1/2 exp (-^) • (1-7) 

Our paper is organized as follows. In Section [2j we recall the replica method and 
the mapping to the 5-Bose gas. To provide some background, we explain the one-point 
generating function for the partition sum and the corresponding Fredholm determinant. 
The extension to two points is discussed and the long-time limit is obtained. The 
technical derivation is carried out in Sections |3] and H] with supporting material in the 
Appendices. 

2. Main results 

2.1. Scale invariance, stationarity 

To have a short hand the Brownian motion average in fOj) is denoted by E^. To 
define it we first introduce the Gaussian average E' 37 with mean and covariance 

E^{x{t)x{t')) = (/3 7 )- 1 min(r, r') , (2.1) 

which corresponds to the Brownian motion starting at 0, and set 

E^ t) (-)=E^(-6(x(t)-x)). (2.2) 

In particular 

,(*(*(<) - x)) = (Pl/2nty^ exp (-^) • (2.3) 

Let us denote the partition function (11.11) by Zp }Jj D(x,t) to indicate explicitly the 
parameter dependence. From the scale invariance of white noise, V(ax, bt) has the 
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same distribution as (a6)~ 1 / 2 V A (x, t), and of the free directed polymer, x(at) has the 
same distribution as a 1//2 x(t), one obtains 

Zf)„A x > V = P^D Z w (/3 3 7j D x, P^D 2 1) . (2.4) 

Thus, it suffices to consider the case when all parameters are equal to one, and from 
now on we adopt the convention 

= 1, 7 = 1, D = 1 , Z(x,t) = Z ltL ,i(x,t) , Eg t) = E ((B , t) . (2.5) 

The free energy ( 11. 5ft has a systematic upward curvature of x 2 /2t, compare with (11.71) . 
However the distribution of F(x,t) — x 2 /2t is independent of the position x. This fact 
can be seen by performing a linear change of variables in the functional integration (11.11) 
defining Z(x, t). One obtains 

(Z(x + a,t))- 1 Z(x + a,t) = (Z{x,t))-% x , t) {exp ( - J &tV(x{t) + ar/t, r)) ) 

= (Z(x,t))- 1 Z(x,t), (2.6) 

where in the second equality we used that the white noise V(x,t) is statistically 
translation invariant in the spatial argument. Hence the distribution of Z(x, t) / (Z(x, t)) 
is independent of x. More general, the stochastic process x i— > Z(x,t)/ (Z(x,t)) is 
stationary in x, and correlation functions of the form 

{fi(Z(x u t)/(Z( Xl , t))) . . . f n (Z(x n , t)/(Z(x n , *)))) (2.7) 

are invariant under a global translation of all the arguments. In particular, the joint 
distribution of Z(xi, t)/ (Z(x\,t)) and Z(x 2 , t)/(Z(x2, £)) depends only on the separation 
x 2 - x 1 . 

For our computation of the two-point generating function it will be of advantage 
to keep the dependence on xi, x 2 separately. The final result will confirm the parabolic 
free energy shift and the dependence on x 2 — x\ only. 



2.2. Replicas 

The n-point correlation function of Z(x, t) can be computed by introducing the replicas 
xi(t), . . . , x n (r), which are simply independent copies of the free directed polymer x{r). 
More specifically, using the explicit form of the generating function for a Gaussian, 

n „i 

(Z(x u t)...Z(x n ,t)} = (l[E {Xj , t) (exp( y - / dry(ar i (r),r)))) 

j= i J o 

n „t n 

= (lI E fe.*))( ex P (s / dr E ^(t) -^(r)))) , (2.8) 
j=i Jo &j=i 

where in the second line the average is over all replicas, the j-th replica starting at 

and ending at Xj at time t. It should be noted that the literal Gaussian average would 

also include the self-interaction term — | Y17=i ~ Xi )- Thus, for the moments of Z, 

the free energy renormalization needed to properly define (11.11) simply corresponds to 

subtract the self-energy. 
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The Feynman-Kac formula implies that the n-point correlation function 
(Z(xi, t) . . . Z(x n , t)) satisfies the imaginary time Schrodinger equation 

- d t (Z(x u t) . . . Z(x n , t)) = H n (Z(x u t) . . . Z(x n , t)) (2.9) 

with the initial condition Z(xj,0) = 8{xj). Here H n is the Lieb-Liniger quantum 
Hamiltonian of n particles on the line with an attractive 5-interaction, 

.. n 1 n 

Hn = -- $>,J 2 - 2 E ~ x i) ( 2 - 10 ) 

[TTl HB] . In this representation the free energy renormalization corresponds to the normal 
ordering of H n . "Solving" (I2.9p . the n-point correlation function is given by 

(Z( Xl ,t)...Z(x n ,t)) = (x u ...,x n \e- tHn \0) . (2.11) 

Here |0) is the state where all particles are at and \xi, . . . ,x n ) is the one where the 
j-th particle is at Xj. Since |0) is symmetric under the exchange of particle labels, one 
can symmetrize in the final state. Thus the propagator e~ tHn is needed only in the 
symmetric sector and the replicas are expressed by the attractive 5-Bose gas on the line. 
As first shown by Mc Guire [18], the ground state energy E {n) of H n is given by 

£oH = -^(n 3 -n) • (2.12) 

The term linear in n translates to the free energy shift t/24, which in fact equals the 
bulk free energy per unit time. For later use, we introduce the parameter 

a = (t/2) 1 / 3 . (2.13) 
In the lowest order approximation, ignoring all excited states, 

(0\e- tHn \0) ^ e~ tEo(n) . (2.14) 
The cubic term of E (n) leads to the decay of the left tail of F(0,t) as 

Prob (F(0, t) - Jjt < u) w exp [-§ (a|u|) 3/2 j (2.15) 

for u — > — oo and t _1 / 3 |w| = O (1) [13 ED], see also [9]. This result confirms that the 
fluctuations of the free energy are of order t 1//3 . In fact, the tail behaviour (12 . 1 5[) agrees 
with the exact tail, see |21|. 



2. 3. One-point generating function of the free energy 

To go beyond (12.141) . one needs the excited states of the attractive 5-Bose gas. They 
can be computed from the Bethe ansatz as has been recently worked out in great detail 
by Dotsenko and Klumov [221 El- I n brackets, we remark that on a ring, the complex 
momenta are solutions of the nonlinear Bethe equations on which little information is 
available. Already to determine the ground state energy requires ingenious computations 
[23J. A corresponding situation has been found for the asymmetric simple exclusion 
process. On a ring, while the ground state and the large deviations for the current have 
been extensively investigated [2U [251 ESI EZ] , the Bethe equations for excited states have 
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been analyzed only partially [2B1 EH GEU EE]- in contrast, for the infinite lattice there 
is a reasonably concise formula for the transition probability with any given number of 
particles [32] . 

Using the complete eigenfunction expansion for H n , a particular generating function 
can be expressed as a Fredholm determinant (T2j [131 E] • More precisely, we define 

G(s; x, t) = (exp (-e^e"^) ) . (2.16) 

Then G is equal to the Fredholm determinant 

G(s-±t-±x 2 ;x,t)=det(l-M). (2.17) 

The operator M does not depend on the position x and depends on time t only through 
the parameter a defined in ( 12.131) . M acts on L 2 (IR) and has the integral kernel 

gait— s 

(u\M\v) = (u\K\v) . (2.18) 

Here K is the Airy operator with integral kernel 

POD 

(u\K\v) = dz Ai(u + z)Ai(z + v) , (2.19) 
Jo 

called Airy kernel. K is related to the Airy Hamiltonian 

H = -(d u ) 2 + u } (2.20) 
as K projects onto all negative eigenstates of H . In particular K* = K and K 2 = K, 



sec 



Appendix B| for details. One easily checks that tr|M| < oo. Hence the Fredholm 
determinant in (I2.17P is well defined. 

The mathematical status of (I2.17P is somewhat tricky. One cannot simply verify 
(I2.17P as the solution of some equation. The derivation relies on the replica method. 
Since log(Z(x,t) n ) ~ n 3 , the moments do not uniquely determine the distribution of 
Z(x,t). To derive (12.171) . one is forced to work with divergent series and has to make 
a reasonable choice for the analytic extension of (Z n ), n 6 N, to the complex plane 
[T2| H3] . However, the generating function G(s;x,t) fixes the distribution of F(x,t), 
which is known by other means [101 E] Thus a posteriori one can verify directly that 
(I2.17P is a valid identity. 

Eq. (I2.17P together with (12.181) establishes that F is of order a. Rescaling s as era 
and taking a — > oo, we note that the right side of (12.161) and the multiplicative prefactor 
in ( I2.18P both converge to a step function. Hence in the long time limit one obtains 

lim Prob (F(x, t) - ±t - ±x 2 < -aa) = det (1 - P a KP a ) = F 2 {a) , (2.21) 

where P a projects on [a, oo). In the long time limit, the free energy fluctuations are 
thus of order t 1//3 . The function F 2 is the celebrated Tracy- Widom distribution function 
[33], which first appeared as the distribution for the maximal eigenvalue of large random 
Hermitian matrices in the Gaussian unitary ensemble (GUE). 

By a more sophisticated argument [32] one also deduces the finite time probability 
density from ( I2.17P with the result 

/oo 
du exp (-e a{s - u) ) g t {u) , (2.22) 
-oo 
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where 

g t (u) = det (1 - P u (B t - P Ai )P u ) - det (1 - P u B t P u ) . (2.23) 
The operators Pm and B t are defined respectively by 

(z\P M \z') = Ai(z)Ai(z') , (2.24) 

and 

roc 1 

(z\B t \z')= I dv- — Ai(z + v) k\{v + z') . (2.25) 

J -co -L e 

There is another model for which the one-point generating function is available [34J: 
one replaces the Brownian motion by a continuous time random walk on the lattice Z 
with forward jumps only. The white noise is correspondingly discretized in the spatial 
direction. The partition function is Zp(N,t) with the polymer starting at and ending 
at (N,t), N > 0. Note that this model has less scale invariance than our case because 
of the lattice. The long time behavior of the free energy remains to be studied. 

For the directed polymer at zero temperature, results are available for a lattice 
discretization in case the random potential has either a one-sided exponential or 
geometric distribution [35] . 



2.4- Two-point generating function and long time limit 

Our novel contribution is the extension of (12.171) to two reference points by using the 
replica method. In analogy to (12.161) . let us define the generating function 

G{s u s 2] x u x 2 , t) = (exp {-e' s ^e- F ^ - e^e"^ 2 '^) ) . (2.26) 

From the mapping of the directed polymer to the £-Bose gas, the generating function G 
can be expanded in a sum over eigenstates of the Lieb-Liniger Hamiltonian f 1 2 . 1 1) . Using 
results from [TJ], under a natural factorization assumption the sum over the eigenstates 
can be written again as a Fredholm determinant. However even then it is difficult to 
extract any useful information from this representation. By a sequence of miraculous 
transformations we arrive at an alternative expression for the Fredholm determinant, 
which turns out to be rather similar to (12.181) in structure and from which the long 
time limit can be read off easily. To distinguish from the true generating function, we 
introduce the sharp superscript " for generating functions with factorization assumption. 
Let us first define the function $ by 

= - (2.27) 

v ' l + e u + e v K ' 

and the operator Q through the kernel 

(u\Q\v) = - si, av - s 2 )(u\e- {2a2) ' llxi - X2lH \v) . (2.28) 

Then 

G*(si Tt x l s 2-2\t- yA, x u x 2 , t) = det (1 - Q ^T^-MHjq _ (Z29) 
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Recall that ( u |e-( 2a2 >~ 1 l' B1 - !Ba l ir |T;) is the propagator of the Airy Hamiltonian and note 
that 

POO 

\ e {2c?)-i\xi-**\H K \ v j = / dze- {2a2) ~ llxi - X2lz Ai(u + z)Ai(z + v) . (2.30) 

Jo 

In particular, e^ 2 " 2 ** X \ X ^- X ^\ H K is a bounded operator. 

From (I2.29P it is obvious that for long times F scales as a ~ t 1//3 and the two-point 
distribution has a non-degenerate limit only if x scales as a 2 ~ t 2 / 3 . So let us substitute 
Si by aa, s 2 by ab, and introduce 

\xi — x 2 \ = 2a 2 y . (2.31) 

Then as t — > oo the left hand side of (12. 29f) converges to the characteristic function 
of the rectangle [— aa, oo) x [— ab, oo), while the operator Q converges to P a e~' y ^ H + 
e -\y\H p b _ p a Q-\v\ H p b . With y, a and b held fixed, we arrive at 

lim Prob [F(xx,t) — — ^x\ > —aa, F(x2,t) — — > —ab) 

t— >oo 

= F 2 {a,b;y), (2.32) 

where the function F 2 (a, b; y) is given in terms of the Fredholm determinant 

F 2 {a, b; y) = det (l - (P a + e'^ H P b e^ H - P a e- m P b e m ) K) . (2.33) 

The function F 2 (a,b;y) is two-point distribution of the Airy process in a form written 
down first in Eq. (5.8) of [21], where the two-point distribution function of the height in 
the polynuclear growth droplet model was studied. For the single step growth model the 
corresponding result was achieved in |36j. Thus apparently the factorization becomes 
exact in the long time limit. 

As shown in [21] the function F 2 (a,b;y) can also be expressed as the Fredholm 
determinant of a 2 x 2 operator kernel, 

^^-(^)(^ K ^J^)}. (,34, 

This form arises naturally when one studies the top most line in Dyson's Brownian 
motion. In the large N limit it converges to the Airy process, which can be viewed as 
the top line of an underlying extended determinantal random field. In particular, the 
n-point distribution is defined most directly through an operator with a n x n matrix 
structure. While we were searching for a corresponding matrix structure, it came as a 
surprise that the replica route apparently prefers the expression (I2.33p . 

From ( 12.341) one can read of properties of F 2 (a,b;y). Obviously it is symmetric in 
a and b. In the limits y — > and y — > oo, the expression of F 2 (a, b; y) simplifies and one 
obtains 

lim F 2 (a, b; y) = F 2 (min(a, b)) and lim F 2 (a, b; y) = F 2 {a)F 2 (b) , (2.35) 

with F 2 of a single argument denoting the Tracy- Widom distribution (12.211) . Using 
(12T33D and [H, K] = one finds 

lim F 2 (a, b;y) = and lim F 2 (a, b; y) = F 2 (a) . (2.36) 

b— >— oo b— >oo 



Two-point function of the free energy for a directed polymer 



9 



Since F 2 (a, b; y) is a distribution function, see (I2.32p . it is an increasing function of a, b. 
It tends to when a, b — > — oo, and to 1 when a, b — > oo. 

3. Replica Summation 

3.1. Two-point generating function 

We start from the expression (12.26P of the generating function G(si, s 2 ; x%, x 2 , t) and 
define 

G\ — G [si — ^jt, $2 — ^t; Xi, x 2 , t) 

= <exp(-e 4/2 V Sl Z(xi,t)-e* /24 e- S2 Z(x 2 ,t))) , (3.1) 

subtracting the linear order of the free energy. Expanding the exponential, G\ is written 
in terms of n-point correlations of the partition function. After decomposing as 

Xi + X 2 X 2 - X! 

x = — g — ' X = — 2 ' ^ ' 

we obtain 

°° ( lWJJV/24 . ^ . 

Gl = 1 + J2 E (ne~ i[(1 ^ i)si+(1+CTi)s2] ) 

iV = l ' <T1,...,0- ]V = ±1 1=1 

x (Z(X + (Tix, t) . . . Z(X + ajvx, t)) . (3.3) 

We use the expression (12.111) for the iV-point correlations in terms of the 5-Bose gas. 
These matrix elements can be expanded as a sum over the orthornormal basis of 
eigenstates, tp r , of the Lieb-Liniger Hamiltonian with the result 

(Z{x 1 ,t)...Z{x N ,t)) = ^2e- tEr (x u ...,x N \^ r )(A\0) , (3-4) 

r 

where E r is the corresponding energy eigenvalue. Combining (13. 3p and (13 .4p . one obtains 
the following expression for the generating function Gi, 

Gi = 1 + £ E e"^ |W(0, • • • , 0) | 2 

AT=1 ~ ' r 

V f TT e-^ft-^)'^^)"^ ^ r(X + ^ • • • X + aNX ^ (3 5) 

^ A 11 / t/ao m ' 1 ' ; 



x 

a\ ,...,<TjV=±l 1=1 



3. ,2. Summation over the eigenstates of the 5-Bose gas 



We have to perform the summation over the eigenstates in the expression (13.51) of the 
generating function G\, for which we follow [H] with some improvements. To ease the 
comparison our notation will be as close as possible to the one in [14] . 

The eigenfunctions of the iV-particle attractive £-Bose gas are labelled by the 
complex wave numbers £ , a — 1, . . . , N. To write them down, we pick positive integers 
n a , a = 1, . . . , M < N such that 

M 

Y^n a = N. (3.6) 



a=l 
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We also introduce the running indices r a = 1, . . . , n a and set 

a 

n («) = J2np, n (0) =0, n^ = N. (3.7) 

/3=1 

Then, for arbitrary vectors q = (qx, . . . , q^) 6 R M and n = (rii, . . . , n M ), one has 

ia = q a - l -{n a + l-2r a ) for a = n^" 1 ) + r a , (3.8) 

with a = 1, . . . , M. According to [14J, Eq. (B.21), the eigenfunction with labels q, n is 
given by 

JV 

= C q$ s % n (p) II (&»(«*) ~ fp(&) + is S n (^a - art)) exp [i £ p(c) a; c . (3.9) 

peV l<a<b<N c=l 

The normalization constant Cq%} is computed in [TJ]. The sum is over the set T 3 of 
all iV-long permutations and sgn(p) is the signature of the permutation p. ipq^n is 
continuous and symmetric in x\, . . . , xn- Each eigenstate of the 5-Bose gas is made 
up of M clusters, where the en's cluster consists of n a bound particles and has center of 
mass momentum q a . 

In the following transformation we invoke a combinatorial identity (we owe the 
proof to P. Di Francesco). 

Lemma 1. Let f(a,b) be arbitrary complex coefficients and let 

d(^,...,^ N )=j2^(p) n &<«)- U) + /m)) ■ ( 3 - iQ ) 

peV l<a<b<N 

Then D equals the Vandermonde determinant, 

D{Z 1 ,...,( N ) = M J] (&-£&)■ (3.H) 

l<a<fe<7V 

Proof. D is a polynomial of total degree N(N — l)/2 with leading coefficient 
N\ rii<a<fe<Af(^a- ~ £&)• Hence, to prove (13.111) one only has to establish that D is 
antisymmetric, since any antisymmetric polynomial is divisible by the Vandermonde 
determinant. 

It suffices to study the interchange of a specific pair, say £i and £2- Let p be some 
permutation such that p(c) = 1, p(d) = 2 and c < d. Then the product in (13.101) 
decomposes into a factor C(p) independent of £1,^2 and a second factor as 

sgn(p)C(p)(6 - £ 2 + f(c, d)) J] (£ p(a) - £1 + f(a, c)) J] (£ p(a) - £ 2 + /(a, d)) 

a<C a < d 

x II & - + 6 )) II & - + /(<*, &)) • (3.12) 

(6 # d) 
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Let p be the permutation with 1 and 2 interchanged. Then C(p) = C(p), sgn(p) = 
— sgn(p), p(c) = 2, p(d) = 1 and the decomposition as in (I3.12p reads 

- sgn(p)C(p)(6 - £i + f(c, d)) J] (e p( a) - 6 + c)) J] - 6 + f{a, d)) 

a<C a < d 

(a ^ c) 

II (6 - e P (6) + f(c, &)) n (6 - C P (6) + /(<*, 6)) • (3.13) 



x 



6 > c 6>a! 

(6 d) 



It is now obvious that the sum of the two terms is antisymmetric in £i and £2, and hence 

D(£i,6,... ) 6v) = -£(k,£i,---,6v)- □ 

We first compute ipq$(0) by taking the limit e — > of x,- = ej. Then sgn(x a — a^) = 
— 1 for all a < 6 and by Lemma CD 

V4#(0) = Cg?J\T! J] (3.14) 

l<a<6<7V 

Secondly we have to evaluate (13.91) at Xj = X + x<jj and to perform the sum over 
all "spin" configurations cr = {<7i, ...,a^}, for which purpose we use that (13. 9 p can be 
written more compactly by using the special structure of the complex wave numbers, 
see [H], Sect. B.2. We introduce the cluster counting function a : [1, N] — > [1, M] 
by 

a(a) = /3 for n (/3_1) < a < n il3) , /3 = 1, M , (3.15) 
and the /3-th cluster by 

fi /3 (p) = {a|a(p(o))=^}. (3.16) 
Then, working out the derivatives in [TJ], Eq. (B.28), 

lf>£S( Xl , Xjy) = CjJ?^ Sgn(p) Yl (?«(p(a)) - ?a(p(6)) + tyfca, ^)) 

pgP 1 < a, f> < JV 

a (p( a )) # <*(?(*>)) 

M , M 

x exp 

a=l cefi a (p) * a=l c,c'eO a (p) 

Here the sum over permutations is understood modulo permutations inside each cluster, 
as indicated by ', and 77 is defined by 

r](x a , x b ) = sgn(x a - x h ) + - ( ^2 s g n ( x « _ x c) ~ ^2 s S n ( x & ~ x <^) ■ ( 3 - 18 ) 

~" c G Q a (p) c 6 n Q (p) 

We spread the particle configuration as 

x a = X + x<j a + ea . (3.19) 

Then 

N 

a i=l 

= E II (9«(PW) - + ^ ■ ( 3 - 2 °) 

(T pgP 1 < a, b < N 

a (p( a )) ^ a(p(6)) 



Af - M 

[i^^ 2 ^"^E l x c-x c /| . (3.17) 
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The phase 4>(<r,p) is given by 

M 

,P) = ^ ( - fOl + S l) n a + - S 2 )m a ((T,p) 

+ iq a (Xn a + xm a ((T,p)) - \\x\{n 2 a - m a ((T,p) 2 )) , (3.21) 



0=1 



where we introduced 

m a (a,p) = a c . (3.22) 

cgO a (p) 

Inserting in ( 13.201) one arrives at 

N 

Yl (n e_ ^ [(1 ^ )Sl+(1+,T!)S2l )^^ X + <?ix + e,...,X + a N x + Ne) 

a i=l 

= C £ ) EE sgn ^ EI (qa( P (a)) - q a ( P (b)) + iv(°a + ea, a b + eb)) 

a peP 1 < a, b < N 

°>(p( a )) # a (p( b )) 

M 

X H I I exp [ - |(si + s 2 )n a + |(si - s 2 )m a (cr,p) 

a=l c£Q a (p) 

+ \q a (Xn a + xm a ((T,p)) - \\x\(n 2 a - m a (cr,p) 2 )] . (3.23) 
Let us shorthand the right hand side of Eq. (I3.23j) as 

5252fi{<r,p)M<r,p). (3.24) 

By Lemma 1 

E/i( CT ^)=^(°) ( 3 - 25 ) 

per 

not depending on cr and, since f 2 (cr,p) depends on cr only through the m Q (cr,p)'s, 
correspondingly 

^/ 2 (cr,p) = c (3.26) 

with c not depending on p. Unfortunately we could not discover any further 
simplification. To proceed anyhow a natural step is to factorize ( 13.231) either with 
respect to p or with respect to cr, both leading in approximation to 

N 

]T ( J]e-^( 1 -^) S1+ ( 1+ ^^])(^(X + cnx + £,... ,X + a w x + Ne)/^{0)) 



U 1=1 

M 

22 II eX P [ - 2"( S 1 + S 2) n a + K S 1 ~~ S 2 )m a ((T 
cr a=l 

+ iq a (Xn a + xm a (cr)) - \\x\{n 2 a - m a (cr) 2 )} , (3.27) 
where 

m a ((T) = ^ °"n(«-l)+r Q • ( 3 - 28 ) 
r a =l 

For small cluster sizes we checked that (13.271) is indeed not a strict equality. 
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4. Two-point generating function 

4-1. Linearization and "spin" summation 

We linearize the terms quadratic in n a and in m a (cr) in the exponential of (I3.27p . so to 
be able to perform the summation over the a^s and the n a 's. For this purpose we use 
the identity 

^au+bv+cuv _ ^cd a d b ^au+bv j\ 

which can be checked by expanding both sides of the equation as a formal power series 
in c, and obtain 

N 

a i=l 
M 

- II e~ x9ld2 exp [ - sx(n a - m a {cr)) - f (n a + m a (<r)) + ig a (Xn Q + OT a (<r))] . 

(4.2) 

Here we introduced the convention 

<9i = <9 S1 , <9 2 = <9 S2 , (4.3) 

which will be used onwards. We note that the exponential inside the product over a 
depends only on the a a with n^ a ~^ < a < n^ a '. Thus, the summation over the a a can 
be performed independently for each a. Recalling ( 13. 2p . we find 

M 

II e~ x9ld2 ex P [ - f K* - m a (tr)) - f (n Q + m a (cr)) + iq a (Xn a + xm Q (cr))] 

cr a=l 

M 

_ J~J ^ix 1 q a -s 1 _|_ e ix 2 i3 Q -S2 j " (4 4) 

a=l 

Fredholm determinant 

We now return to the generating function Gi of (I3.5p . denoting by its approximation 
under the factorization assumption (I3.27p . The eigenfunctions ipq$ are normalized in 
such a way that they form an orthonormal basis in the symmetric subspace of L 2 (M N ). 
With this normalization, and using (13. 14ft . 

l^ } (0,---,0)| 2 = iV!det( . * a) ) , (4-5) 

\5V n i + n k) +l{Qj - Qk)J jjk =l,...,M 

see Eqs. (B.58) and (34) of [Hj. From Eq. (B.29) of [H], the energy E ( q % ] of the 
eigenstate if)q$ is 
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Finally, the properly normalized sum over the eigenstates is given by 

E^sn(/jE)w... 

see Eqs. (B.53) and (B.60) of [H]. Combining all, we obtain for G\ the following 



expression 

00 ^ M 



«l-^EBrn(r*5:)-fi 



1 



M 

t3 , 



i=i 



We observe that ( 14.81) can be rewritten as a Fredholm determinant (see Appendix A for 
a few basic facts about Fredholm determinants). Indeed, if one introduces the kernel R 

as 

R(q, n- q', n') = w , a T •/ * — , (4-9) 

27r ±(n + n'J + i(q - q') 

then the generating function G\ is given by 

G\ = det(l + R) . (4.10) 

More explicitly, it holds 

Gf = 1 + / d ?i • • • d( ?M J] det ( R (lj> n i'i Ik, n h )) j>k=lj ^ M . (4.11) 

M=l ' ni,...,n Af =l 

We perform the summation over the n/s and the integration over the q/s. For this 

purpose the integrated version for the denominator in (14. 9p is used, 

1 /*°° i 

r^— -=/ dze- z ^ n+n ' )+iiq -^ . (4.12) 

|(n + ri) + i(q - q') J 

The operator R can then be written as a product of two operators, R = R1R2, 



00 



R(q,n;q',n')= / dz -Ri(g, n; z)R 2 (z; q',n') , (4.13) 

—00 

with 

#i(g, n; z) = t {z>0} e-^V^V^ 1 ^ - e -§* e -V/2 + e i* 2? - S2 ) p ? ( 414 ) 



and 



it! 2 (z; g',n') = 1 { , >0} ^ e^'V* . (4.15) 



Using det(l + R1R2) = det(l + R2R1), the generating function G 1 } becomes equal to a 
Fredholm determinant of the new operator N, 

G\ = det(l + N) , (4.16) 
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where N = R2R1 with kernel 

N(z, z' 



/oo ^ 
dg^i2 2 (^; ?) n) J R 1 (g,n;y) . (4.17) 

00 71=1 

The variables qj and n,, which were previously the variables corresponding to the 
definition of the Fredholm determinant, are now inside the kernel N. 



In Appendix C , the summation over n and the integration over q in (14.171) is 
performed explicitly. Most of the steps are rather similar to the computations done 
in [13], [12j [T3] in case of the one-point generating function. Note that on face value the 
sum over n is badly divergent. In terms of the parameter a = (t/2) 1 / 3 (equal to 2 2 / 3 A 
in the notation of [H]) and of the function $ defined in equation (I2.27p . the kernel of 
iV equals 

N(z, z') = -a~ l N(a- l z, a~ l z) , (4.18) 

with 

/oo 
-00 

$(mt - St, au - s 2 )Ai(z + u)Ai(u + z) . (4.19) 

The generating function G\ is now given by 

G\ = det(l - N) . (4.20) 

We will simplify the kernel and express it in terms of the Airy Hamiltonian H and 
of the Airy operator K. 

4-3. Subtraction of the parabolic profile 

The following transformations are guided to have the shift by x 2 /2t manifestly visible in 
G\. Using the definition (13.21) of X and x, the expression (I4.19P of the kernel iV rewrites 

as 

N(z,z') = l {z ^ >0} e- xdld2 e^ idl - d2) ^- d ^ 

/oo 
du e -|( 9 i+*)( 9 -^')$( au _ Su au _ S2 )Ai(z + u)k\[u + z') . (4.21) 
-00 

We note that d\ + 82 in (I4.2ip can replaced by — a -1 <9 u , where the derivative d u acts 
only on $(cra — s\, au — s 2 ) and not on the product of Airy functions. One can then 
make d u to act only on the product of Airy functions by integrating by parts. Thereby 

N(z, z') = t^^e-^e^ 91 - 82 ^- 9 ^ 

du §(au - si, au - s 2 ) e - (2a2) " lx ^ (9z -^' ) Ai(^ + u)Ai(u + /) . (4.22) 



Since the Airy function is a solution of the differential equation Ai"(n) = uAi(w), we 
have 

d u (d z - d z ,)Ai(z + u)Ai{u + z') = {z- z')Ai{z + u)Ai(u + z') (4.23) 
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and hence 

e a[d u (d z -d z ,)-(z-z')] Ai ( z + U ^ M ( U + ^ = + U ^ M ( U + ^ _ ^ 

The commutator of <9 U (<92 — d z i) and — z') is equal to 2<9 U and it commutes with both 
d u (d z — d z >) and (z — z'). Thus the Baker-Campbell-Hausdorff formula terminates and 

e -adu(d z -d z ,) Ai ( z + M )Ai( M + = e - a(z - y) e a29 "Ai(z + u)Ai(u + z') . (4.25) 

We use this property in the expression of the kernel N, and integrate by parts to make 
d u act again on $(cra — si, au — s 2 ) with the result 

/oo 
d M e-( 2Q2 ^ lx ( 2 - z Ui(^ + M ) 
-oo 

x Ai(« + z')e-^~ lx2d ^(au - Sl ,au- s 2 ) . (4.26) 
Using the commutation relation 

e a(d z -d z ,) e -b(z-z') _ e -2ab e -b(z-z') e a(d z -d z ,) ^ (4.27) 

one obtains 

/oo 
duM(z + u) 
-oo 

x Ai(« + ^)e-( 2a3 )" x ^ 9l - 92 )e-( 4a4 )" x2 ^$(a M - Sl ,au- s 2 ) . (4.28) 
Since e ad acts as a shift operator, we arrive at 

N(z, z') = 1 {2 , 2 , >0} e -( 2 « 2 )-^( z - z ')e"^ 92 e^( 9l - a2 )^-^)e-( 4Q3 ^ la; ^ 9l+92 ) 

/oo 
duA\(z + u)M(u + z 1 ) <E» (au - si - ^x\,au - s 2 - Yt x l) ■ (4.29) 
-oo 

The factor exp ( — (2a 2 )~ 1 X(z — z')) can be eliminated by a similarity transformation 
of the kernel, it does not contribute to the Fredholm determinant. We define the shifted 
generating function 

G 2 (u, v; x, t) = G\(u — ^x\,v — ^.xl; x±,x 2 , t) . (4.30) 
and find that G\ can be written as 

G\ — det(l - L) , (4.31) 
where the kernel L is given by 

L(z,z') = i {ZiZ , >0} e-^ 92 e^^- 92 )^-^)e-( 4a3 ^ lK2 ^ +92 ) 

/oo 
du M(z + u)Ai(u + z') <E» (au - si, au - s 2 ) . (4.32) 
-oo 

In this form, the kernel depends only on x (and t), as to be expected from the discussion 
in the introduction. 
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4-4- Rewriting of the kernel L in terms of the Airy Hamiltonian 

The next step is to eliminate the operator did 2 from the expression ( 14. 32ft for the kernel 
L, for which purpose one writes 

/oo 
-oo 

x $(cm - si, au - s 2 )Ai(z + u)Ai(u + z') . (4.33) 

In this expression, we can replace (d\ + d2) 2 by a~ 2 dP, where d u acts only on 
$(cm — s\,au — s 2 ) and not on the product of Airy functions. Then, integrating by 
parts, d u acts on the product of Airy functions instead. Thereby 



L(z,z') = 1 {Z , Z , >0} / due^^e-^^^e^ 81 -^ 9 '- 9 ^ 
J —00 

x - si, aw - s 2 )e- (2a2) " lx9 «Ai(2 + u)Ai(u + z') . (4.34) 

But d u acting on A\(z + u)Ai(u + z') is the same as d z + d z >. Thus, one can replace 
exp( — (2a 2 )~ 1 xd^) by exp( — (2a 2 )~ 1 xd u (d z + d z /)). After integrating by parts again, 
we replace exp((2a 2 ) -1 :r<9 u (c^ + d z >)) by exp( — -§-(di + d 2 )(d z + d z >)) and obtain 

/oo 
-oo 

x $(cra - si, aw - s 2 )Ai(z + w)Ai(it + z') . (4.35) 

We now introduce a dummy integration to achieve that au — s±, au — s 2 , Ai(z + u) and 
Ai(w + z') do not all depend on the same variable u. It holds 

/oo 
dudv6(u - t ,) e i(«-f^) e 5«-S*«.)e-(*» 8 )- li »'(*+*) 
-00 

x - si, aw - s 2 )Ai(z + v)Ai(u + z') . (4.36) 

Note that 

[(dl - 2a~ 1 d 1 d z/ ) + a' 2 (z' + H u )] $(a« - s u av - s 2 )Ai(z + v)Ai(u + z') = , (4.37) 

where H u = — <9 2 + u denotes the Airy Hamiltonian acting on the variable u. The 
commutation relation 

\(d\ - 2a~ 1 d 1 d z ,), a-\z' + H u )\ = -2a^d 1 (4.38) 

together with the Baker- Campbell-Hausdorff formula implies 

e ¥ 9 l-^ d z') $(q, m _ Sl) av _ s 2 )Ai(z + v)Ai(u + z) 

A corresponding relation is obtained by interchanging the roles of z, Si, u and of z', s 2 , 
v. Using both the kernel L becomes 

/oo 
du dv 6{U - v ) e -i^ 2 )-^{z + z') e -{2o?)-^{H u+ H v ) 
-00 

x $(au - si, av - s 2 )Ai(z + v)Ai(u + 2;') . (4.40) 
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A final integration by parts over u and v yields 

/oo 
dudv (u\e- a ~ 2xH \v) 
-oo 

x $(cra - si, aw - s 2 )e- {2a2) ' lx{z+z,) Ai(z + v)Ai(u + z') . (4.41) 

This last expression is an operator product of the form ABA. Hence using the cyclicity 
of the determinant one arrives at 

G\ = det(l — L) — det(l - L) , (4.42) 

where 

(u\L\v) = (u\e- a ~ 2xH \v}$(au - s u av - s 2 )(u\e a ~ 2xH K\v) . (4.43) 
We conclude that L = Qe a 2 \ X \ H K, as to be shown. 

5. Finite time probability density function 

The two-point free energy fluctuations can be written as 

nxj,t) = -^t+±x} + $i(t), J = 1,2, (5.1) 

with random amplitudes £,j(t). In the factorization approximation, we know already that 
£j(t) = O (t 1 / 3 )- The joint distribution depends only on \xi — x 2 \, and nondegenerate 
correlations occur for a separation of order t 2//3 . Following the procedure in [12], we 
would like to extract the underlying pdf from the generating function G\ Let us denote 
by pl{wi, w 2 ) the approximate joint pdf of £i(t) and £, 2 (t), and write it as the convolution 
with two independent Gumbel densities Fq u : 

/oo 

dvj. dv 2 Fg u (wi - v 1 )Fq u (w 2 - v 2 )g t {v x ,v 2 ) , (5.2) 
-oo 

where 

F Gu (w) = exp(-e- w ) , (5.3) 

compare with the one-point distribution (12.221) . Eq. (15.21) defines the yet to be 
determined function g t , which is normalized to 1 by construction. From numerical 
solutions in the one-point case, we know that gt of f !2.23f) is in general not everywhere 
positive. This implies that g t of (15. 2p will not be a pdf, in general. 
The generating function G\ reads 



G\= \ dw 1 dw 2 p t (w 1 ,w 2 )exp(-e- Sl - Wl -e- S2 - W2 ) . (5.4) 

J — oo 

Inserting (15.21) and performing the integration over w\ and w 2 yields 

G\ = / dvi dv 2 qt(vi,v 2 ) . (5.5) 

2 / 2utyi, z) eVl + e _ Sl e „ 2 + e _ S2 v ) 
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We analytically continue on both sides from e _Sj to — e aj — iajE, e > 0, j = 1, 2. From 
the identity 

lim >^ = 2in5(v — a) , (5.6) 

cr=±l 

the left side of f !5.5j) multiplied by o"iO"2 and summed over o~i, a-i = ±1 yields in the limit 
e -> 

- 47r 2 ^(ai,a 2 ) • (5.7) 

On the right side of (I5.5p . using G2 = det(l — L) with L given by (14.411) . we obtain the 
sum of four Fredholm determinants with operators L aiU2 , o~j — ±1. To compute these 
kernels, we have to take the limit e — > of 

duidu2h(ui,u 2 )- —, 7, (5.8) 

with h general at this stage. Using 

v(-)+unr6(y), (5.9) 



y - iae \y 

one obtains 



/oo 
dudv (u\e- a ~ 2xH \v)e- {2a2) ' lx{z+z,) Ai(z + u)Ai(v + z') 
-00 



_gQ«l+ai gOU2+02 



x (_ e «n+ai _ e aU2+a2 )5(l - e QUl+ai - e a " 2+a2 )) , (5.10) 

and 

gt(a 1 ,a 2 ) = -(2tt) -2 ^ aio- 2 det(l - L CTl>(T2 ) . (5.11) 

CTl,(T2=±l 

In the one-point case, inserting the corresponding analytic continuation in (I2.18P 
yields a one-dimensional projection. This further simplifies the expression for the pdf, 
compare with (I2.23p . For the two-point case, no further simplification seems to be 
available. 



6. Conclusions 

Recently the probability distribution of the free energy of the point-to-point continuum 
directed polymer has been computed exactly, using the approximation through the 
weakly asymmetric simple exclusion process [101 E] • This result could be reproduced 
by using replicas and the complete eigenfunction expansion of the propagator of the 
attractive 5-Bose gas on the line [121 EE] . 

In our contribution we studied the joint pdf of F(xi, t), F(x 2 , t). In this case 
the approximation through the weakly asymmetric exclusion process, while still valid, 
is no longer computable and we have to rely on the replica method, which yields a 
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particular generating function. Invoking a specific factorization, the result is expressed 
as a Fredholm determinant. In fact, the corresponding operator has a structure rather 
similar to the case of one point, compare (12.1 7ft and (12. 29 h . 

Equipped with this information, we established the large time limit of the pdf 
yielding a result in agreement with lattice directed polymers at zero temperature. In 
1 + 1 dimensions all models with a short range disorder potential are expected to flow to 
the zero temperature fixed point. Our result further supports this claim. 

The free energy of the point-to-point continuum directed polymer is isomorphic 
to the solution of the KPZ equation with sharp wedge initial data. Thus we have 
automatically determined the joint pdf of the heights h(xi,t), h(x2,t) of the KPZ 
equation for large times, in particular the height-height correlation function. This 
function has been measured recently for droplet growth in a thin film of turbulent liquid 
crystal [37]. The experimental curve agrees very well with the theoretical prediction in 
the limit t — >• oo. 
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Appendix A. Fredholm determinants 

Let us first consider the case of a finite n x n matrix A. The Taylor expansion of the 
determinant of 1 + zA is given by the von Koch formula [38] in terms of the minors of 
A. It holds 

n m n 

det(l + zA) = J2^J E det (%Oi<^< m - (A.l) 

m=0 ' ii,...,i m =l 

If the matrix A is now replaced by an integral operator A with kernel A(u,v ) = (tt|A|u), 
von Koch formula (lA.lj) formally rewrites as 

det(l + zA) = ^2 — r / d«i • • • du m det (A(uj, «i0)i<j,fc<m • (A. 2) 

m=0 J 

Of course, du could mean a more general summation procedure. In particular, it could 
refer to summation over some discrete index and integration over K. flA.2j) can be used 
as the definition of the Fredholm determinant det(l + zA). In order for this definition 
to make sense, the operator A is required to be trace-class, we refer to [39J for details. 
If so, the logarithm of the Fredholm determinant is given by 

logdet(l + zA) = trlog(l + zA) 

= E / dui ■ ■ ■ dUm A ( u u u 2 )A(u2,u 3 ) . . . A(u m , «i) • (A.3) 

T77 / 
m=l J 
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Another useful identity for Fredholm determinants is the cycle property. If A and B are 
Hilbert-Schmidt operators [i.e. ti AA* < oo and trBB* < oo), then 



This property allows to exchange the roles of integrations which are inside the defining 
kernel of AB with the integration corresponding to the Fredholm determinant. We 
emphasize that the two kernels AB and BA do not necessarily act on the same space. 

Numerical evaluations of Fredholm determinants can be performed by discretizing 
the integrals in ( 1A.2I) . Using the von Koch formula ( lA.lj) . the evaluation of a Fredholm 
determinant is thereby reduced to the computation of the determinant of a finite matrix. 
We refer to [38l 140] for an illuminating discussion and precise error estimates. 

Appendix B. Airy operator and Airy Hamiltonian 

We recall the definition of the Airy Hamiltonian H and of the Airy operator K. We 
work in the space of complex- valued square integrable functions L 2 (1R) with the scalar 
product 



If necessary, we write H as H u to indicate the variable on which the Airy operator is 
acting. The Airy function is the solution of the differential equation 



det(l + AB) = det(l + 5A) . 



(A.4) 




(B.l) 



The Airy Hamiltonian H is defined by 
H = -(d u f + u. 



(B.2) 



AS" {u) = uM(u) 
such that Ai(u) — > as u — > oo. Setting 
(j) z {u) = Ai(u - z) , 



(B.3) 



(B.4) 



one notes that 4> z satisfies the eigenvalue equation 



H<fi z = z(f) z . 



(B.5) 



In addition, 




(B.6) 



In Dirac notation this completeness relation reads 




(B.7) 



Hence the Airy Hamiltonian has the spectral representation 




(B.8) 
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The projection to all negative energy states defines the Airy operator 

K= f dz . (B.9) 

J — oo 

In particular, one has K = K*, K 2 = K, and obviously [K,H] = 0. In position 
representation, the Airy kernel writes 

/ , r „ v f°° , » / w ., x Ai(u)Ai'(u) - Ai'(u)Ai(u) ,„ , 

(uK v) = / dzAtiu + z)AUz + v) = — ^ — — — — — . B.10 

Jo u-v 

We also define the projection onto the spatial interval [a, oo) by 

POO 

P a = du\u)(u\ . (B.ll) 



The operator P a KP a is trace class for all a > — oo. F 2 (a) = det(l — P a KP a ) is by 
definition the Tracy- Widom distribution [33] corresponding to the Gaussian Unitary 
Ensemble of random matrices. 

Appendix C. Integration over q and summation over n in the kernel N 

We start from the explicit expression 

N(z,z') = t {z , z , >0} ^fv«(-V« 3 / 2 V^ 

J -co n=1 

X [ - e -|(*+*') e -t<z72 ( e ixig--i + e ix 2 g-s 2 ^ jn _ (q-q 

In order to perform the summation over n and the integration over q, we insert the 
classical relation 

POO 

e tn3/24 = / duAi ( u )eW 1/3nu . (C.2) 



Besides the Airy function there are infinitely many other functions which satisfy (1C.2j) . 
One concrete example would be 

Ai{u) + sin(vTM)e- u2 / 2 . (C.3) 

Our choice is determined by being the only one which provides the correct one-point 
result. Using the binomial theorem to expand the term of power n, one obtains 



N(Z.Z') = 11;,,.. „; ^ '^f2J2^b¥ q[Z 



x f ^_ e (t/8) 1 /3 Ue _i(, + ,') e _ t g2/ 2 j- e -x9 1 9 2e - Sl fc- S2 (n-fc) _ ^ ^ 

We introduce the parameter a = (t/2) 1 ' 3 and perform the change of variable u — > 
u + 2 2 / 3 a 2 q 2 + 2~ 1 / 3 a~ 1 (z + z') in the integral. This results in 

N(z, z') = 1 {Z<Z , >0} ^ jr Ai (« + 22/3 « 2 ? 2 + ^cr\z + z')) 

J -°° ZlX n=l fe=0 
X e' 9 ^ - «'+a:ifc+£C2(n— fc)] | J ^_^ n e 2 ~ 2/3aun Q- x 9id2 e -sik-s 2 (n-k) ^ 
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Using the relation 

r S Ai(ag2 + b)eicq 

= 2" 1 / 3 a- 1 / 2 Ai (2- 2 / 3 (b + a-^c) ) Ai {2^ (b - aT^c) ) , (C.6) 
see [4T] , the integration over g can be performed. One obtains 

00 n =l k=0 ^ ' 

x 2~ 2/3 cr 1 Ai (2~ 2/3 w + oT x z + (2a) -1 [xiJfe + x 2 {n - k)}) 

x Ai (2- 2/3 u + a~V - {2a)~ 1 [x 1 k + x 2 {n - k)]) . (C.7) 

We change variables as u — > z — >■ az, and use the relation 

f(z + a)=exp[ad z ]f(z) (C.8) 

to move n and k out of the Airy functions. In preparation for the summation over n 
and k, one finds 

„oo 00 n , v 

aiV(a;z, c^) = l {z , z ' >0} / du ^^^g-^ft^^-^tn-t) 

n=l fc=0 ^ ' 

x e (2a)- 1 ^ifc+^(n-fc))(9,-9,0 Ai ( u + 2 ) Ai ( M + . (Q.9) 

Noticing the identity 

e (2a)- 1 (x 1 k+x 2 (n~k))(d z -d z ,) e ~s 1 k-s 2 (n-k) _ ^-(2a)- 1 (x 2 d 2 +x 1 d 1 )(d z -d z ,) e -s 1 k-s 2 (n-k) ^q^ 

the summations over k and n finally yield the expression (I4.18p . (14.191) for the kernel N 
with <3> is defined as in Eq. (|2.27j) . 
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